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where 1^9 is a vector function on R'' belonging to some Banach space V, X isa stochastic 
process and the integral is some version of a stochastic integral defined via regular- 
ization. A stochastic current is a continuous version of this map, seen as a random 
element of the topological dual of V. We give sufficient conditions for the current to 
live in some Sobolev space of distributions and we provide elements to conjecture that 
those are also necessary. Next we verify the sufficient conditions when the process X is 
a d-dimensional fractional Brownian motion (fBm); we identify regularity in Sobolev 
spaces for fBm with Hurst index H G (1/4,1). Next we provide some results about 
general Sobolev regularity of Brownian currents. Finally we discuss applications to a 
model of random vortex filaments in turbulent fiuids. 
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Abstract 



We study the pathwise regularity of the map 




1 Introduction 



We consider stochastic integrals, loosely speaking of the form 



Iiv)= [ {ipiXt),dXt) 



(1) 



Jo 



where (Xt) is a Wiener process or a fractional Brownian motion with Hurst parameter H 
in a certain range. We are interested in the pathwise continuity properties with respect 
to ip: we would like to establish that the random generalized field ip ^ I (ip) has a version 
that is a.s. continuous in ip in certain topologies. In the language of geometric measure 
theory, such a property means that the stochastic integral defines pathwise a current, with 
the regularity specified by the topologies that we have found. 

This problem is motivated by the study of fluido dynamical models. In in the study 
of the energy of a vortex filament naturally appear some stochastic double integral related 
to Wiener process 



where f{x) = (x) where {x) is the kernel of the pseudo-differential operator (1 — 
A)^" (precise definitions will be given in section [2^ . / is therefore a continuous singular 
function at zero. 

The difficulty there comes from the appearance of anticipating integrands and from 
the singularity of / at zero. [1] gives sense to this integral in some Stratonovich sense. 
Moreover that paper explores the connection with self-intersection local time considered 
for instance by J.-F. Le Gall in [lOj . 

The work [T6j considers a similar double integral in the case of fractional Brownian 
motion with Hurst index H > ^ using Malliavin-Skorohod anticipating calculus. 

A natural approach is to interpret previous double integral as a symmetric (or even- 
tually) forward integral in the framework of stochastic calculus via regularization, see [20j 
for a survey. We recall that when X is a semimartingale, forward (respectively symmet- 
ric) integral (p{X)d~ X (resp. ip{X)d° X) coincides with the corresponding Ito (resp. 
Stratonovich) integral. The double stochastic integral considered by [1] coincides in fact 
with the symmetric integral introduced here. So ^ can be interpreted as 



In this paper, X will be a fractional Brownian motion with Hurst index H > 1/4 but 
a complete study of the existence of integrals Q will be not yet performed here because 
of heavy technicalities. We will only essentially consider their regularized versions. Now, 
those double integrals are naturally in correspondence with currents related to / defined in 
([1]) . The investigation of those currents is strictly related to "pathwise stochastic calculus" 
in the spirit of rough paths theory by T. Lyons and coauthors, [151 114j . Here we aim at 
exploring the "pathwise character" of stochastic integrals via regularization. A first step 
in this direction was done in [12j, where the authors showed that forward integrals of the 




(2) 




(3) 
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type Jq ip{X)d~ X , when X is a one-dimensional semimartingale or a fractional Brownian 
motion with Hurst index H > can be regarded as a.s. uniform approximations of their 
regularization I~{^p) (see section 3), instead of the usual convergence in probability. 

This analysis of currents related to stochastic integrals was started in using an 
approach based on spectral analysis. The approach presented here is not based on Fourier 
transform and contains new general ideas with respect to [9]. Informally speaking, it is 
based on the formula 

r {if (Xt) , dXt) = [ /(I - A^if (x) , r (x - Xt) dXt) dx. (4) 

Jo jRd \ Jo I 

This formula decouples and X and replaces the problem of the pathwise dependence of 
/ on the infinite dimensional parameter Lp with the problem of the pathwise dependence 
of Jq^ Ka [x — Xt) dXt on the finite dimensional parameter x G M*^. Another form of 
decoupling is also one of the ingredients of the Fourier approach of [9j but the novelty here 
is that we can take better advantages from the properties of the underlying process (like, 
for example, the existence of a density). Moreover formula ([4]) produces at least two new 
results. 

First, we can treat in an essentially optimal way the case of fractional Brownian motion, 
making use of its Gaussian properties. For H > 1/2 results in this direction can be 
extracted from the estimates proved in [T6] again by spectral analysis. However, with the 
present approach we may treat the case H G (1/4, 1/2) as well. 

Second, in the case of the Brownian motion, we may work with functions ip in the 
Sobolev spaces of Banach type Hp, with p > 1, instead of only the Hilbert topologies H2 
considered in [9], with the great advantage that it is sufficient to ask less differentiability 
on if (any a > 1 suffices), at the price of a larger p (depending on a and the space 
dimension). In this way we may cover, for instance, the class G C^'^ treated in |15j . 
see also |13j : the approach here is entirely different and does not rely on rough paths, see 
Remark [371 

Finally, we apply these ideas to random vortex filaments. In the case of the fractional 
Brownian motion we prove new results about the finiteness of the kinetic energy of the 
filaments (such a property is expected to be linked to the regularity of the pathwise 
current). In the case of the Brownian motion optimal conditions for a finite energy were 
already proved in [31[5], while a sufficient condition when H > 1/2 has been found in |16j . 
The results of the present work provide new regularity properties of the random filaments, 
especially for the parameter range H £ (1/4, 1/2). 

2 Generalities 

2.1 Stochastic currents 

Let (Xf) be a stochastic process such that Xq = a.s. on a probability space {i^,J^,P) 
with values in M"^. Let T > 0. Let y be a Banach space of vector fields : M"' ^ M"' and 
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D C V he a dense subset. Assume that a stochastic integral 

I{^)= f {^{Xt),dXt) 
Jo 

is weh defined, in a suitable sense (Ito, etc.), for every ip £ D. Our first aim is to define 
it for every ip £ V. In addition, we would like to prove that it has a pathwise redefinition 
according to the following: 

Definition 1 The family of r.v. has a pathwise redefinition on V if there 

exists a measurable mapping : 0, ^ V' such that for every (p £ D 

I (ip) (io) = (io)) {ip) for P-a.e. io £ Q. (5) 

Then, if we succeed in our objective, 

• for every ip £V we consider the r.v. a; i-^ (u)) (y?) as the definition of the stochastic 
integral / (y?) (now extended to the class ip £ V ) 

• for P-a.e. u; S 17, we consider the linear continuous mapping ip (^(cj)) (ip) as a 
pathwise redefinition of stochastic integral on V. 

Formally, the candidate for ^ is the expression 

e(x)= r 5{x-Xt)dXt. 
Jo 

where 5 is here the d-dimensional Dirac measure. Indeed, always formally, 

i{^)= I {C{x),^ix))dx= r (( [ 6{x-Xt)^ix)dx) ,dXt) 

jRd Jo \ VJRd / / 

= / {ip{Xt),dXt) = Ii^). 
Jo 

We remark that this viewpoint is inspired by the theory of currents; with other methods 
(spectral ones) it was developed in [S]. 



2.2 Decoupling by duality 

As we said in the introduction, our approach is based on a proper rigorous version of 
formula One way to interpret it by the following duality argument, that we describe 
only at a formal level. 

Let W be another Banach space and A : V ^ W he an isomorphism. Proceeding 
formally as above we have 

/ {ipiXt),dXt)^, = {^,Ovy' = (^"'Me)yv" 
Jo 

= (a^, (a-i)*A 

\ ^ ^ /w,w' 
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(notice that A ^ : W ^ V, (^A : V' ^ W'). Our aim essentially amounts to prove 
that (A-^)*? : n^W is a well defined random variable. 

This reformulation becomes useful if the spaces W, W' are easier to handle than V, 
V, and the operator (A~^)* has a kernel K{x,y) as an operator in function spaces: 



((A~i)*/) {x)= j K{x,y)f{y)dy. 



In such a case, formally 



((A-i)*e) ix) = I K{x,y) 6{y-Xt)dXt^dy 

= £ (^j K{x,y)6{y- Xt) dy^ dXt 

= j K{x,Xt)dXt. 
Jo 

Below we make a rigorous version of this representation by choosing V = Hp (R"^) , W = 
LP {R'^), A = (1 - A)t, K = (x - y) (notations are given in the next section). 

2.3 Rigorous setting 



Denote by S (M j the space of rapidly decreasing infinitely differentiable vector fields 
99 : R'^ — > M'^, by S' (W^^ its dual (the space of tempered distributional fields) and by T 
the Fourier transform 

(T^) {£) = (2^)-'^/2 f e-'^'-"^^ip{x)dx, eeR'^ 



which is an isomorphism in both S (M'^) and S' (R'^). Let J-~^ denote the inverse Fourier 
transform. For every s E R, let A^ : S' (R'') S' (R*^) be the pseudo-differential operator 
defined as 

s_ 



We shall also denote it by (1 — A) 2 . 

Let iW^) , with p > 1 and s G M, be the Sobolev space of vector fields f €^ S' 
such that 

M\ms ■= \(1 - A)iip(xWdx < 00, 



see |21j . sec. 2.3.3, where the definition chosen here for brevity is given as a characteri- 
zation. From the very definitions of A^ and Hp (R'') , the operator Ag is an isomorphism 
from Hp (R*^) onto (R*^) (the Lebesgue space of p-integrable vector fields). 
Another fact often used in the paper is that the dual space [Hp (R*^))' is Hp,^ 
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see [21] . section 2.6.1. Moreover, being an isomorphism from (M'^) onto (K'^), 



its dual operator A* is an isomorphism from (M j onto H~i' 

It is known that negative fractional powers of a positive selfadjoint operator A in a 
Hilbert space H, such that —A generates the semigroup T{t), have the representation 



1 r°° 

A-^ = —- r-^T{t)dt, 

r (a) Jo 



where F is the standard Gamma function, see [17J, formula (6.9). Taking A = (1 — A) in 
the Hilbert space H = (R°') , we have 



(T it) ip) (x) = (47rt)-'^/2 / e-^-V(2/)d2/ 

and thus, for a > 0, 

((1 - A)-» (x) = ^^^^ r [ e'-^-'^{y)dydt 
r(a) Jq J^d 

oo 



1 /"OO , I 



r (a) (4^)'^/^ Jo 



2 e 4t 



In fact this formula can be proved more elementarily from the definition of (1 — A) 
and the formula 

1 



oo 







then taking A = and the Fourier transform of the Gaussian density. This fact implies 

that the operator (1 — A)~", which originally is an isomorphism between H2 ° and 
(M'^), considered by restriction as a bounded linear operator in has a kernel 

((l-A)-»(x)= / K^{x-y)ipiy)dy (6) 

given by 

Ka{x) = TTT / r 26 4t * — . 

F (a) (47r)'^/2 7q t 

The following estimates are not optimized as far as the exponential decay is concerned; 
we just state a version sufficient for our purposes. The proof of the two lemmas before are 
in Appendix A. 

Lemma 2 There exists positive constants Ca^d; Ca^d such that: 
1 ) For < a < I , we have 



Ka (x) = p{x) (7) 



where Ca^e ^'^1 < p (x) < Ca^e 
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2) For a > we have 

_l£li _ M 

Ca^e 4 < Ka (x) < Ca4e 8 (8) 

for two positive constants c'^^, C'^^. Moreover 

K^{x) = K^O) - p'{x)\x\-^+^^ > (9) 

with < Ka{0) < oo and where Ca^rfe"^'^'' < p' {x) < Ca,de~^~^ ; 

3) Finally, when a = d/2 we have 

i^a(x) <C,,d loglxle-"^" N (10) 

where is another positive constant. 

Remark 3 In particular p and p' are bounded. 

In the applications we will need also some control on AKa{x) which is provided by 
the next lemma. 

Lemma 4 It holds that —AKa{x) = Ka{x) — Ka-i{x). Then, when a < d/2 + I, 

|-Ai^„(x)| </(x)|x|2"-^-2 (11) 

where p" is positive, bounded above, locally bounded away from zero below and depends on 
a. 

2.4 Regularity of stochastic currents 

With these notations and preliminaries in mind, we may state a first rigorous variant of 
formula Given a continuous stochastic process {Xt)^^Q on {Q,J^,P) with values in 
W^, given e > 0, let {DeXt)^^^ be any one of the following discrete derivatives: 

Xt+e — Xt Xt — Xt-e ^t+e — Xts 

see 
where we understand that Xt^^ = for t < e. The following integral 



Ie{^) = I {^{Xt),D,Xt)dt 

Jo 



is well defined P -cl.S. 3S Si classical Riemann integral, at least for every continuous vector 
field if. 

Lemma 5 Given a,e > 0, with probability one the function t Ka/2 ~ ^t) De^t is 
integrable for a.e. x G M*^, the function 



r?,(x):= / K^/2{x-Xt)D,Xtdt 
Jo 



is in L} (M'^) and for any (p G S (R*^) we have 

r {if {Xt) , DeXt) dt= f /(I - A)"/V {x) , r {x - Xt) D.Xtdt) dx. (12) 
Jo Jr'^ \ Jo I 

Proof. Notice that (1 - A)"/^ g S (M"^) and, by ([6]), 



/ K^l2^.-x)\{\-^rl\ 



{x)dx. 



Thus 



^ {^{Xt),D,Xt)dt= r If K^i^iXt-x) (1-A)-/V {x)dx,D,Xt)dt. 
Jo \Jr'' 



Denote by (x) the function equal to {x) for x / 0, infinite for x = 0. Suppose 
for a moment that 



P 



{Xt — x) dxdt < 



oo 



1. 



(13) 



Then the integrabihty properties stated in the lemma will hold. Since 



sup 

{t,x)e[0,T]xI 



1-A)-/V {x),D,X, 



< oo U.S., 



using Fubini theorem we will get (fT^ (notice that i-^t — x) = (x — ^t)). 

Thus we have only to prove p^ . Since Ka/2 is positive, we may apply again Fubini 
theorem and analyze J^iJ^d Ka/2{Xt — x)dx)dt. But we have, for every y 



Ka/2 {y - X) dx 



Ka/2 (x) dx. 



This quantity is finite for every a > 0, from the estimates of lemma [2j The proof is now 
complete. ■ 

Below we need a criterion to decide when tj^ (defined in the previous lemma) belongs to 
CW^) . It is thus useful to introduce the following condition which ensures the existence 



of the representation given in Lemma [8] below. 
Condition 6 {Xt — Xs) dtds < oo P-a.s. 



Remark 7 /q' Ka {Xt - Xs) dtds < oo, implies that the function {t, s) ^ Xt — Xs 
is different from zero except possibly on a zero measure set of [0,T]^, and that the well- 
defined function {t, s) ^ Ka {Xt — Xs) is Lebesgue integrable on [0, T]^. From now on the 
notation will simply be replaced by K^- 
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Lemma 8 Under Condition\^we have the following double integral representation for the 
norm of r]^ defined in Lemma\^ 



i2(K.) = Ka {Xt - Xs) {DeXt, D,Xs) dtds. 



Proof. We have 

|2 



||r/e||^2(Kd) = / {r}e{x) ,r]e{x))dx 

^ LXl ^"/2 - ^ K^/2{x-Xs)D,X,ds^dx 

f [ [ K^/2{x-Xt)K^/2{x-X,){DeXuDeX,)dxdtds 
Jo Jo Jm.'^ 

if we can apply Fubini Theorem. Then it is sufficient to use the property 

Ka/2 {x - y) {x - z) dx = {v - z) . 



Since the process X is continuous, both (t, s) ^ Df,XsDf,Xt is a continuous two-parameter 
process which on [0, T]^ is a.s. bounded. Then a sufficient condition to apply Fubini 
Theorem is 



ill {x - Xt) {x - Xs) dxdtds < oo. 

Jo Jo Jr'^ 

Since the integrand is positive, it is equal to 

[ [ ([ K^/2{x-Xt)K^/2{x-Xs)dx)dtds= [ [ K^{Xt-X,)dtds. 
Jo Jo kJm.'' / Jo Jo 

Invoking condition [6] we can conclude the proof. ■ 

The double integral representation of the norm of ??e will play a major role in the 
following, so we introduce the notation 

Za^e := [ [ Ka (Xt - Xs) {DeXt, D^Xs) dtds. 

Jo Jo 

Lemma 9 Assume Condition\^ holds for any a>a. Then the function a i— > Z^^e >0 is 
decreasing for a >a. 

Proof. Denote rja^eix) = Jq Ka/2{x — Xt)Di,Xtdt making explicit the dependence on 
a. It is not difficult to prove that, if a < a < /? we have rjj^^^ = (1 — A)("~^)/2T7a^£. 
Then Z^,^ = \\r]p^e\\\2^^d) = ||(1 - A)("-^)/2r?c,,e||^2(Rd) < \\riaA\\2(Rd) = ^o.,e, since being 
a — /? < , the operator (1 — A)^""^)/^ has a norm bounded by one. ■ 

We will assume below the following condition on the convergence of the regularized 
integrals. 
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Condition 10 For every if G {^'^) , h (v?) converges in probability to some r.v., de- 
noted by I {if). 



Under Condition [TUl the mapping ip ^ I (ip) is a priori defined only on (M j with 
values in the set ($7) of random variables. Its extension to 99 G H2 (M"^) is a result of 
the next theorem. 

Theorem 11 Assume Conditions and \l(A and the a priori bound 

cT f-T 



sup E 

£6(0,1) 



f [ Ka{Xt-Xs){DeXt,DeXs)dtds 

Jo Jo 



< 00. 



Then: 



i) the mappings (p G (M'^) ^ ((/?),/ ((/?) G (Q.) take values in and 
extend (uniquely) to linear continuous mappings from {R'^) to (Q). Moreover, for 
every (p G H2 {R'^) , le — > I {'■p) in probability and in L'^~^ [VL) for every (5 > 0. 

ii) In addition, there exist random elements : 17 (W^^ (in fact belonging 
to L^(il; //^"(M'^)) j that constitute pathwise redefinitions of and I over the functions 



, in the sense of Definition\^ 



Proof. Step 1 (mean square results). By the assumptions and the previous lemma 
we have sup£g(o,i) ^ 



I l|2 



< 00. From (HI]), for V? G {W^) we have 

/.M = ((i-Ar/2. 



L2 



Therefore, always for 93 G C| 

\Ie (<^)| < ||%|Il2 

E 



;i-A) 



a/2 



L2 



I l|2 



<cM\l, 



Immediately we have 1^ (v?) G {^) for every ip G (M*^) , and the mapping ip 1— > 
Is (92) extends (uniquely by density) to a linear continuous mapping from H2 {W^) to 

L2 (0). 

Given <p> G Cq°(M'^), since lei'p) — > in probability, uniform integrability arguments 
and E'[|/e((/?)p] < C4||(/^||^a(]gd); yield hiip) I{ip) in L^~^(il) for every 6 > 0; moreover, 

Tf.„\ r- T'2ln\ I?rir/',oM2l ^ ll,„l|2 

HS 



it is not difficult to deduce I{ip) G L^(Q) and £'[|I(93)p] < C'Q||v9||^aj^jgd^. As before, this 
implies that the mapping ip ^ I [ip) extends uniquely to a linear continuous mapping 
from H2 (M'^) to (0). Now, with these extensions, it is not difficult to show that 
Is (ip) — > / (ip) in L^^^ (Q) for every 6 > also for every for ip G H2 



Step 2 (pathwise results). We still have to construct and ^. Recalling that ??e G 
L2 (0;L2(R'^)), is simply defined as [(1- A)"/2]*r/e, element of (j^jiJ^"" (M'^)), where 
A* denotes the dual of an operator A. 
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To this end, recall that (1 — A)"^^ is an isomorphism between H2 (M*^) and L^l 
and thus the dual operator [(1 — A)"/^]* is an isomorphism between the dual spaces 
and i?2~" {^'^) (we identify L'^{W^) with its dual). 

The family {r]e} is bounded in (il; L^(IR'^)), hence there exist a sequence r/e^ weakly 
convergent to some r/ in (il; L'^iW^)) : E {ri^„,Y) j^2^d'^ E {r], Y)^2(^d-) for every Y € 

L2(0;L2(]R'^)). We set ^ := [(1 - A)°/2]*r/, random element of H^"" (M"^). We shall see 
that this definition does not depend on the sequence e„. We have to prove that for every 
if G Cq° (M*^), / (99) (1^) = (lo)) (if) for P-a.s. a; € Jl. Equivalently we have to prove that 
for every ip e {M.'^) 



M=r?(a;) ((1 - Af/^) 



for P-a.s. w E 



(14) 



We already know that leif) = ((1 — ^)"'^^V5) ^e)L2(]Rd) for every e > 0. Choose Y above 



of the form Y = F {I - A)"/^ with generic F £ L'^ (n). Given ip £ C^ 
that 



^ (M'*) , we know 



F(r?,„,(l-A) 



a/2 



L2 



F/r?, (1-A)"/^ 



L2 



for every F £ (0). Hence 



F(??,(l- A 



,a/2 



L2 



but we also know that ((^) ^ / (93) in ^ (il) for every (5 > 0. We get 



E[FI{lp)] = E 



F/r?, (1-A)"/^ 



L2 



at least for every bounded random variable F, hence (I14p holds true. This also implies 
that the definition of 5, does not depend on the sequence e„. The proof is complete. ■ 

To state a possible converse of Th. [TT] it is useful to introduce a weaker version of 
Def.E 

Definition 12 Let A £ be such that P{A) > 0. We say that {I{ip)}y:,QV has a pathwise 
redefinition on A if there exists a measurable mapping ^ : A ^ V' such that for every 
if G {W^) equation ([^ holds true for P-a.e. lv £ A. 

Theorem 13 Assume Conditions and \10\. and the a priori bound 



E 



Jo Jo 



Xs){D,Xt,D,Xs)dtds 



< 00 



for every e > 0. If there exists A£ T with P{A) > such that {I{p)}^ev ^o,^ ^ pathwise 
redefinition on A, then 

limsup I I Ka {Xt - Xs) {DsXt, D^Xs) dtds < 00 P-a.s. lo G A. 
£^0 Jo Jo 
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Proof. The first part of step 1 of the previous proof is still valid (except for the 
uniformity in e of the constants). Thus in particular i?[||%||^2(iRd)] < cxd, Is{v) = ((1 — 
A)°/V,r/,), \I,{ip)\ < C^JifWl^^^.y Ie{v) e L^{n) for every ^ G Co°°(M^) and the 

mapping Lp ^ ((/?) extends to a linear continuous mapping from H2 (M'^) to (17). 
We know, by Condition [10] and the existence of ^, that for every S (M°') 

{l-/\r/\,r,,{u:)) ^ (e(^))(^) forP-a.e. e A 

One can find a countable set T> C Cg° (M*^) with the following two properties: i) (1 — A)"*^^ T) 
is dense in L^(M'^) and ii) for P-a.e. uj & A 

l{l-Ar/\,7j,{u:)) ^ (e(^))(<^) for every G P. (15) 

Let us prove the claim by contradiction. Assume there is A' £ J^, A' C A, P{A') > 0, 
such that 

limsup / / Ka (Xt — Xg) {D^Xt, D^Xg) dtds = 00 for every uj G A' . 
e^o Jo Jo 

By lemma [SI limsup ||r/e (w)||^2(iRd) = 00 for every a; G A'. 

e—>0 

Consequently, there is a subset A" of A' with P{A") > such that (|15p holds true for 

every to G A" . Thus, given lo G A" , there is an infinitesimal sequence {e„ (i^)} such that 
II 1 1 2 

lim„_oo Pe„(a;) {^)\\ i^2^^d-^ = oo but at the same time 

hm ((1 - A)"/V,^.„H H)...o^ = 

for every G P. This is impossible because of the density of (1 — A)"''^ P in The 
proof is complete. ■ 

Definition 14 PFe say that a.s. a stochastic current I does not belong to V' if there is 
no A £ P{A) > such that {I[ip)}^p^v has a pathwise redefinition on A. 

Corollary 15 Under the conditions of Theorem \13l if 

limsup [ I Ka (Xt - Xs) {D^Xt, D^Xs) dtds = +00 P-a.s. 

e-^O Jo Jo 
the current I does not belong to H2 
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3 Application to the fractional Brownian motion 



We recall here that a fractional Brownian motion (fBm) B = {Bt) with Hurst index 
H G (0, 1), is a Gaussian mean-zero real process whose covariance function is given by 

Cov{B,, Bt) = + If ^ -\s- {s,t) G R^. 

This process has been widely studied: for some recent developments, we point to [3] as 
a relevant monograph. For instance we recall that when H = ^, B is a. classical Wiener 
process. Its trajectories are Holder continuous with respect to any parameter j < H. 
Recall that 

E{\Bt-Bs\') = \t-sf''. 

For our stochastic integral redefinition, we have chosen the framework of stochastic calculus 
via regularization; for survey about the topic and recent developments, see [20] . 

Let X = [X^, ...,X'^) be a d-dimensional fractional Brownian motion with Hurst index 
H G (0,1), i.e. an M'^-valued process whose components X^,...,X'^ are real independent 
fractional Brownian motions. In this section we study the regularity of the current gen- 
erated by X using the results of the previous section and in particular Theorem [TT] which 
gives sufficient conditions for regularity in the Hilbert spaces (^^'^)- 

We will consider the symmetric and forward integrals, respectively defined as the limit 
in probability, as e — > 0, of 



and 



:=£ <^iXt),^^^±^^^>dt 



I-{ip) := I < ipiXt), ^^l±L_Jk > dt. 
'0 £ 



Whenever they exist we will write 

r{ip) = / < ^{Xt),d''Xt >, and I-{ip) = / < ip{Xt),d-Xt > . 
JO Jo 

In [1] the authors show that the symmetric integral exists for fBm with any H > 1/4 
(the proof about the d = 1 case but it extends without problem to higher dimensions). For 
H > 1/6 necessary and sufficient conditions for the existence of the symmetric integral 
are given in |ll]. 

As far as the forward integral is concerned, it is known that, in dimension 1, it does 
not exist at least for some (very simple) functions 93. Existence of the forward integral in 
dimension one is guaranteed if and only if > 1/2 [IT9t I20|. Observe that, when H < ^ 

the forward integral Bd^B does not exist since B is not of finite quadratic variation 
process. 

When H > 1/2 and for d > 1 the forward integral (equal to the Young integral in the 
case H > 1/2), is equal to the symmetric integral minus the covariation [ip{X), X]/2. This 
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exists if X has all its mutual covariations and the i-th component of that bracket gives 

d 



f{X),X% = fY.djf{X,)d[X\X^l 



see again [19]. 



Before entering into details concerning stochastic integration, we state an important 
preliminary result. 

Proposition 16 If a > max(0, d/2 - 1/{2H)) then E f^|^ Ka{Xt - Xs)dsdt^ < +oo. 
Therefore Condition holds. 

Proof. Since the corresponding random variable is non-negative, previous expectation 
equals 

/ / E{K^{Xt-Xs))dsdt = 2 f [ E{KjXt-X,))dsdt. (16) 
Jo Jo Jo Jo 

Next we use the bound on given in Lemma [2) 

• Suppose first that d/2 — 1/{2H) < a < |. Then the previous expression is bounded 
by const Jq^ Jq -Ed — Xs\'^"~'^)dsdt. For any 7 > —d, the scaling property of 
fractional Brownian motion gives -E[|Xt — Xs\'^] = -£^[1 A^|'''](i — s)'^^ , where is a 
standard d-dimensional Gaussian random variable. Then the right member of (jl6p 
equals 

const S[|iV|2"-'^ /" [ \t- s\'^'^^-^^dsdt (17) 
Jo Jo 

which is finite if 2aH — dH > —1 and a > 0. Therefore when a > max(d/2 — 
1/{2H),0). 

• Suppose now that ^ < a. In this case Ka < const ()16p is trivially bounded. 

• Finally, in the case a = d/2 we have Ka{x) < clog \x\ and using again the scaling, 
it is easy to prove the boundedness of ([16]). 



The Proposition above will allow us to verify the condition required by Lemma [8] of 
previous section. 

3.1 Symmetric integral 

Let D^Xt = {Xt+s - Xt^e)/2e and let / : M'^ - {0} ^ M+ such that 

/ I fiXt - Xs)dtds <oo P- a.s. (18) 
Jo Jo 
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and denote ^ ^ 

Zsif) ■■= [ [ f{Xt- Xs) {D^Xt, DlX,)dsdt. 
Jo Jo 

In this section we will study the r.v. Z^^e = Zi;{Ka) in order to obtain necessary and 
sufficient conditions for the regularity of the symmetric current I{ip) based on fBm. 
The following lemma is proved in Appendix A. 

Lemma 17 1. We have the following estimates: 

\Cov{D°,XlD°^Xi)\ < const\t - sp^'^; (19) 

and 

I Cov{DlXlXl - = I Cov{D^^XlXi - Xi)\ < const\t - s\'^"-^ (20) 
uniformly in e > 0. 

2. If s ^ t, one has 

lim I Cov{Dlxl,Xl - Xl)\ = 2H\t - s\'^"-^. 

The main results of this section are contained in the next two theorems. Let an = 
d/2-l + 1/{2H). 

Theorem 18 For any H > 1/4 and any a > an we have sup^ EZa,e < c». 

Remark 19 We observe that d/2 — 1/2H < an so that for the range value for a in 
Theorem [TSi Condition is verified. 

Theorem 20 For any H > 1/4 and any a < an we have liminf^ £'(Za^e) = +oo. 

Remark 21 The statement of Theorem <ind Lemma below allow us to formulate 
the following conjecture that we have not been able to prove: for if > 1/4 and a < an we 
should have 

liminf Z^e = +oo, a.s. 

For {I {'■p)} ^p(^v to have a pathwise redefinition on some A £ J^,P{A) > a necessary 
condition is that limsup£_^o ^a,£ < oo on A (Th. \T3\) . 

If this conjecture were true we could establish that a.s. I{^p) does not belong to ff2^"(M'^) 
when a < an- 

Before proving the Theorems we deduce Sobolev regularity of fBm with any Hurst 
parameter between (1/4, 1). 

Corollary 22 (Regularity of symmetric currents) The symmetric integral of a frac- 
tional Brownian motion with Hurst parameter H > 1/4 admits a pathwise redefinition on 
the space ii^"(]R'^) for any a > ajj. 
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Proof. By Theorem 1181 we know that E[Za^^] is uniformly bounded in e when a > an- 
Since the regularized integrals Is{v) converges in probabihty as e — > for any H > 1/4, 
Condition [TU] holds. So we can apply Theorem II II to obtain a pathwise current with values 
in H^^iM.'^) for any a> an- ■ 

In the following proof we will use a basic result about Gaussian random variables 
recalled here: 

Lemma 23 ("Wick's theorem") Let Z = (Z£)i<£<7v be a mean-zero Gaussian random 
vector and f £ C^(K^ -jM.), then we have 

N 

E[ZJ{Z)] = Cov{Ze, Z,)E[Vjf{Z)]. (21) 
j=i 

Proof. The conclusion follows easily taking /(zi, ztv) = 6xp(z ^^-^ i^Zj) for any t = 
(ii, ...,tN) G K^. In fact, in that case E{f(Z)) is provided by the characteristic function. 
Therefore one has 

/ tTt'^ 

E{ex.p{it ■ Z)) = exp 



where t' stands for the transposition and F is the covariance matrix of Z. Differentiating 
the previous expression with respect to ti provides the result (I2ip for the particular case 
of /. The general result follows by usual density arguments. ■ 

Proof of Th. 1181 Using Lemma [23l independence and equal distribution of different 
coordinates we have 

EZ^^e =E r r (Xt - X,) V Cov{D^,Xi, D^.Xl) dtds 
Jo Jo ~^ 

+ E [ [ AKc, (Xt - Xs) Cov{D^X^ , X^ - X} ) Cov(L>°X] , X} - X} ) dtds 
Jo Jo 

f-T j-T 

= E / KaiXt- X,)y2Cov{D'^^XlD'^^Xl)dtds 
Jo Jo j 

+ E / {-A)K^ {Xt - Xs) \Coy{D^,X},X} - 
Jo Jo 



:l)\'^dtds 



(22) 

where we used the fact that 

Coy{D',Xlxl - Xl) = -Cov{Dlx],Xl - X]) = 1 (|t - . + ep^^ - |t - . - ep^) 
which can be verified by a straightforward computation. 
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Consider first the case H > 1/2. Assume a < d/2 and note that when H > 1/2 we 
have an < d/2. By lemma [T71 we have 



rT rT 

\2H-2 



EZa,e < const E [ I Ka{Xt- Xs) \t- s 
Jo Jo 

+ const / / \{-A)Ka{Xt- Xs)\\t - s 

Jo Jo 



4H-2 



dtds 



to JO 

and then, using lemma [2] 1) and 4), we get 

EZr,, < const E [ f \X,-XJ^''-'^\t-s\^^~^dtds 



■'OL,£ ^ const XL/ I I 

/o Jo 

+ const / / \Xt- Xs\'^°'''^'^\t- s\'^^-'^ dtds 
Jo Jo 

const / / 



Jo 



(23) 



+ const r r \t - s|(2"~'^+2)/^-2 ^^^^ 

^0 Jo 



which are uniformly bounded in e if (2a + 2 — d)H > 1, i.e. when a > an as required. 
We have established the uniform bound when a G {an, d/2) but now recall that EZ^^e is 
a decreasing function of a so that this bound extends to all a > an- 

Let us now consider the case H < 1/2 (so that now an > d/2) and assume a > d/2. 
Rewrite Za,e as 

Za,e = Z,{K) + ZeiK^m (24) 

where ha{x) = Ka{x) — i^a(O). Note that < Ka{Q) < oo when a > d/2. Moreover 

f-T rT 



Z,{K^{Q))= [ [ K^mD'iXs^DlXt) = K^{<d) 
Jo Jo 



DlXsds 



D^Xtdt = / '~' dt = / Xtdt 



and 

lo " Jo £ Jt 

so that by the continuity of the process X we have the limit 

limZ,{K^{0))=K^mXT\^ 

exists almost surely and Zg^KaiO)) is uniformly in (actually in all L^). 

So it remains to consider Zi,{ha)- For we have the estimate Q, provided < 
2q! — d < 2, so that we obtain an upper bound similar to eq. (I23p and the same condition 
on a follows. Note that a must satisfy 

d 

d/2 < an < a < - + l 
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so that we must require H > 1/4. Again by monotonicity of a ^ Z^^^ we have uniform 
boundedness for any a > an- 
m 

Proof of Th. 1201 We will perform a decomposition of Zq,^,. as follows. Write 

where ^ ^ 

Aa,e = [ [ K^iXt- Xs) V CoviD^^Xf, Dlxi) dtds 
Jo Jo 



Ba,e= [ [ i-A)Ka{Xt - Xs)\Cov{D^,X^,Xl - Xl)\^ dtds 
Jo Jo 

and Qa,e is the remainder. Note that, by comparing this decomposition with eq. (j22p . we 
have 

EZa^e = EAa^e + EBa,e 

SO that EQa^p = 0. This is a kind of Wick product decomposition, but not quite, since 
the terms A, B are not constants, but still random variables. 
A useful remark is that A^^e ^ since we can write 

Ao,,e = E [ r {Xt - Xs) {DlXt, Z)°X,) dtds 
Jo Jo 



tribution of X and where E denotes expectation with respect to this auxiliary fBm. So 
we have the formula A^^^ = E\\fiir\\'^ where fi^{x) = Jq K^/2ix ~ Xt)D'^Xtdt which shows 



where we introduced an auxiliary independent d-dimensional fBm X with the same dis- 
tribution of X 
we have the fo: 
that Aa,e > 0. 

Next using the equality —AKa{x) = Ka-i{x) — Ka{x) we rewrite Ba,e as Ba): — Ba}e 
where ^ ^ 

5^=/ / K^-iiXt-Xs)\CoviD^,X^,X^-Xl)\^dtds 
Jo Jo 

4S=/ / K^iXt-Xs)\CoviD^,Xl,Xl-Xl)\^dtds. 
Jo Jo 

(2) 

Let us show first that E\B£ \ is uniformly bounded in e when a > a// — 1 and H > 1/4. 
Indeed when a < d/2 by computations similar to those of Th. [18] we have that E\Bs \ 
is uniformly bounded if a > an — 1 and a > 0. On the other hand, when a > d/2, the 
kernel is bounded, so Lemma [T71 1) allows to write 

I < const / / \Coy{D^X^,X^ - dtds < const / [ \t- s\^"-^ dtds < C 
Jo Jo Jo Jo 

uniformly in e provided H > 1/4. 

Moreover below we will show the following. 
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Lemma 24 If H > 1/4 and a < an then liminis^o Ba^ = +00 a.s. 

Then if we admit the result of previous lemma we can conclude with the use of Fatou 
lemma. In fact, for any a £ (an — 1, C(h) and for some positive constant c we have 

liminf > limini E{Bae) > ^^(liminf S^,^)) - sup E{B^^l) 

e—*0 ' e— »0 ' e^O ' £ ' 

> ^(liminf _ c = +00. 

Moreover observe that this is enough since, from LemmaO we have that Za,e is a decreas- 
ing function of a so the result will hold for any a < an- ■ 



Proof of Lemma | 

By Fatou lemma we have 



liminf B^^) > 2 / [ {Xt - X^) liminf |Cov(D°X/, - X})]^ dtds 

But, when t ^ s, 



liminf |Cov(Z)°X/,Xi^ _ = iH'^\t - s 



AH-2 



. Now assume 

a<d/2 + l. (25) 

then by Lemma[2]there exist a small constant r > such that > C|xp°^'^^^l5(o,r)(^)- 

This allows us to bound from below as follows 

lim inf 5 W > const / / Ko,-i (Xj - X,) |t - s\^"-^ dtds 

> const / / \Xt-Xs\'^'''^'^^^lB(o,r){Xt-Xs)\t-s\'^^-^dtds. 
Jo Jo 

Since the paths of fBm are Holder continuous with parameter strictly smaller than H, for 
any 'j < H there exists a random constant Cx,^ such that 

\Xt-X,\<CxMi-s\\ t,s£[0,T] 

By choosing a random time S > small enough such that sup^ .,g[o,5] l-^t — Xs\ < r we 
have 

rS rt 

liminf S^i) > const / \Xt - X,\'^°'-'^-'^\t - sj^^'^ dtds 
-0 ' s t (26) 

> const C^-^"+2 /■ f 

Jo Jo 

where 5 = {d + 2 — 2a){H — 7) is a arbitrarily small positive constant since 'j < H can be 
chosen arbitrarily near to H and 

d + 2-2a>4-4>0 
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when H > 1/4. Then when a < an we can choose 6 smah enough to make the double 
integral in eq. (p6|) diverge. Summing up we must have a < mm{aH, d/2 + 1) and H > 
1/4. and H > 1/4. But when H > 1/4 we have oh < d/2 — 1 so that sufficient conditions 
are a < an and H > 1/4. This observation concludes the proof. ■ 

3.2 The Forward integral 

Let D-Xt = {Xt+e - Xt)/e, take / : R'^ - {0} ^ M+ satisfying ([18]) and denote 

Zrif) = r r f i^t - Xs) {D-Xt, D-Xs)dsdt. 
Jo Jo 

We can state similar theorems as in previous subsection. 

Theorem 25 For any H > 1/2 and any a > an we have sup^ EZ^^e < oo. 

Theorem 26 For H > 1/2 and a < an we have liminf^ £'(ZQ,^e) = +oo. 

The proofs follow the same line as the corresponding theorems about symmetric inte- 
grals. For this the following lemma will be crucial. 

Lemma 27 1. We have the following estimates: 

I Cov{D-Xl D-Xi)\ < const\t - s\'^^~^ (27) 

and 

I Cov{D-Xi,Xi -Xi)\ = \ Cov{D-Xi,Xi - < const\t - s|^^"^ (28) 
2. If s ^ t, one has 



lim I Cov{D-Xl,Xl -Xl)\ = 2H\t 



e^O 



S 



2H-1 



Proof of Th. 1261 Again the proof is similar to the one of Th. [20] where D~ is replaced 
by D^, Lemma [27] is used instead of lemma [T7] In particular, according to Lemma [271 if 
s one has 

liminf \Cow[D-Xl,X} - X])p = AH'^\t - s\'^^-'^ 

e^O 

when H > ^. ■ 

In particular we have the following Corollary. 

Corollary 28 (Regularity of forward currents) The forward integral of a fractional 
Brownian motion with Hurst parameter H > 1/2 admits a pathwise redefinition on the 
space ff2^"(M'^) for any a > an- 

Proof. By Theorem [25] we know that EZ^^s is uniformly bounded in e when a > an- 
Since the regularized integrals converges in probability as e ^ for any H > 1/2, 

Condition 1101 holds and we can apply Theorem [11] to obtain a pathwise current with values 
in //2^"(M'*) for any a > oh- ■ 
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4 Brownian regularity in Hp, J9 7^ 2 

In this section we restrict ourselves to the case when X is a d-dimensional classical Brow- 
nian motion, that we denote by W. The key ingredient is the following lemma. 

Lemma 29 // the dimension d >2 and the real numbers a > 1 and p' > 1 satisfy 

{d-a + l)p' <d 

then 



E 



exp (— e \x — Wt 



p'/2 



\2d-2a 



'-^dt 



dx < oo 



\x- Wt\ 
for every e > 0. 

We shall prove below this lemma. Let us first describe its consequences. 
From the bounds on K^/2j see [21 we have 

' rT \ P'/21 



J^^E Kl/,ix-Wt)dt^ 



dx < oo 



and thus for a.e. x G M'^ we have PiJ^ K'^^^{x — Wt)dt < oo) = 1 which implies that the 
Ito integral 



T] (x) := 



r ix-Wt)dWt 
Jo 



is well defined, for a.e. x G M.'^, as a limit in probability of 



K, 



a/2 



[x - Wt) dt, 



see for instance |20) . These approximation integrals are measurable in the pair 
hence they are measurable in a; as a mapping with values in the space of random variables 
with the metric of convergence in probability, and this way one can see that the limit object 
r/ (x) is measurable in the pair {x,u:). From Burkhoder-Davies-Gundy (BDG) inequality 
we have E \r]{x)\^ dx < oo and thus r] G L'p (il x M'^) and 



P [uj e Q : X ^ r] {x,uj) E 



1. 



To minimize the subtleties related to a direct use of rj (x), we introduce a regularization. 
Let T (t) be the semigroup on (R'^) generated by (A — 1) and we set 



K 



(S) 



%{x) := {T{5)K^/,) (x) = (4vr5) 



-d/2 



\x-y\^ ^ 

e 4* ''K^/^ {y) dy. 
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We have xfj^ £ S (M^). Set 



T?^'^) {x) 



Jo ' 



which is obviously weh defined for every x and has a measurable version in the pair 
{x,uj). It is not difficult to justify that r]^^'> is square integrable in {x,uj) and that 77^''^ = 
T {5/2) r/(^/2) hence r?^'^) G S {W^) with probability one. 
We have the following regularized version of (jH). 



Lemma 30 For 5 > 

t-T 







{{T{5)^){Wt),dWt)= I Ul-/:^r'\{x)J'\x))dx. 



(29) 



Proof. Given a vector field (f> S (M'^) and a continuous exponentially decreasing 
function ip on M'^, we have the Fubini type identity 



{Wt — x) (f){x)dx, dWt 

lo \Jr'' 

with probability one. We omit the details of the proof. We have 
T{5)^ = T {5) (1 - A)-"/2(l - A)"/V 



and thus 



{{T{5) 99) {Wt),dWt\ 







x)dx, / V (W^t - a;) dWt 



(1 - A)"/Vl 



a/2 



K^l (■ - x) 



a/2 



Here we can apply the Fubini rule because T (6) tp €z S (M"^) for 5 > and (1 — A)"/^c/? G 
S (W^) . This implies ([2^]) and completes the proof. ■ 



(30) 



Lemma 31 For d > 2, a > 1, p' > 1, such that {d — a + 1) p' < d we have 



sup / E 

5>0 



r/(^) (x) 



dx < 00. 



Proof. Let us restrict the argument to the most difficult case < a < d where Ka/2 
has a singularity at zero. From BDG inequality we have 



/ 



E 



T]'^^^ {X) 



dx < Cr, 



E 



dt 



p'/2- 



dx. 



From the definition of K^^^j^ in terms of and estimate ([7]) we get the inequality 



K^;],{x)<C^,,{A7:5)-'"^ 



e 4* \y\ e & dy. 
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Let us show that this impHes 



e 8 



(31) 



for a new constant C^^d, uniformly in 5 E (0,1). The proof of this result for |x| > 1 is 
rather easy, so let us only deal with \x\ < 1. Write x = re with |e| = 1 and change variable 
y = rz in the integral, to get 



K[%{x)<r"-''Ca,d{'^7rS)-''''r 



-d/2 d 



2i |2 
r \e — z\ ^ ^ 

43 \z\ dz 



\a—d 



(see the definition of the semigroup T{t)). It is now easy to see that T (t) 

bounded above by a constant, uniformly in t > 0. This proves ([3T]) . 

Having this estimate, it is sufficient to apply lemma [221 The proof is complete. 
We can now prove the main result of this section. 



(e) is 



rp 

Theorem 32 The ltd integral Jq {(f{Wt) ,dWt) has a pathwise redefinition on the space 
V = H^ (M*^) for every dimension d and real numbers a > 1 and p > 1 satisfying 

d 



p > 



a — 1 



In particular, in any dimension d, given e > 0, for every p > ^ the integral {(f (Wt) , dWt 
has a pathwise redefinition on the space H^^'^ (M'^) . 



Proof. Step 1. In the case d = 1 we have 
theorem (see |21j, section 2.8.1, remark 2). Thus 



C (M) by Sobolev embedding 



^ {Wt) dWt = $ (Wt) - $ (0) - J if' {Wt 
Jo ^ Jo 



dt 



where = (p. This implies the result. We restrict now to the case d> 2. 

Step 2. We pass to the limit in ()29p . Let us treat the left-hand-side. With easy 
manipulations we see that 

(T {5) ip) {x) = (27r)-'^/2 e-^ [ (x - zV26) dz 

Hence, splitting the integral in a sufficiently large ball and the complementary, since ip € 
S (M'^) , we see that T {6) f ^ 'p uniformly over all R'^ as 5 ^ 0. Thus J^^ {{T {6) 99) {Wt) , dWt) 
easily converges to {cp {Wt) , dWt), in mean square. 

Given the value of p in the statement of the theorem, under the assumption a > 1 the 
inequality p > is equivalent to {d — a + l)p' < d, where 1/p+l/p' = 1, so the previous 
lemma applies. From (j30p there is a sequence 5„ ^ and an element rj^^^ G L^' (17 x M'^) 
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such that 77'^''") 



(p ^ S (W^^ , we thus have, in the hmit as 

{^{Wt),dWt 



weakly in ($7 x M'^), when n +00. From (j29p . for a given 



E 



X 







X 

n ■ 
X 



00, 



;i_ A)°/V(x),r/(°) {x))dx 



for every bounded r.v. X and thus 



/" {^{Wt),dWt)= [ ((1-A)"/V(x),7?(°)(x))dx 

JM'' ^ ' 



with probabihty one. 

Step 3. Therefore, given S (M'^), with probabihty one we have 



{piWt),dWt: 



< 



(1 - A)"/V 



(0) 



V 

(0) 



LP 



The proof is complete. ■ 

Remark 33 The same result is true for the stopped Brownian motion 

= WtArn, TR = inf {t > : \Wt\ > R} . 

with given R > 0. The statement is that the ltd integral (99 (M^/^) , dW^^ has a pathwise 
redefinition on the space (M*^) under the same conditions on d,a,p as in the theorem. 
The proof is the same ( even easier, since in the proof of lemma (23 we do not have to care 
of the exponential term). 



Remark 34 The same result is true in the Sobolev-Slobodeckij spaces {R'^) defined 

rT 



JHy, section 2.3. The statement is that the ltd integral Jq {ip{Wt) ,dWt) has a pathwise 
redefinition on the space Wp (M°') under the same conditions on d,a,p as in the theorem. 
Indeed, given a triple d,a,p as in the theorem, let a' < a be such that also the triple 
d,a',p satisfies the assumption of the theorem. Then fQ{(f{Wt),dWt) has a pathwise 
redefinition on H^' (M*^) ; by definition of pathwise redefinition, we see that this implies 

that Jq {(f{Wt) ,dWt) has a pathwise redefinition on the space Wp (M*^), because we have 
the continuous embedding 

( R'^] C H° 



see ]21^ . remark 4 of section 2.3.3. The same result is of course true for the ltd integral 



We can now elaborate the previous results in the direction of the Holder topology. 
Given e e (0, 1), denote by (7^+^ (M'^) the space of all continuously differentiable functions 
/ on R^ such that 



^.i+e = sup (|/(x)| + \Df{x)\) + sup 



\Df{x)-Dfiy)\ 



< 00, 
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see [21], section 2.7. Endowed with the norm the space C^~^^ (M'^) is a Banach 

space. 

Theorem 35 In any dimension d, for every e S (0, 1) the ltd integral {(f (Wt) , dWt) 



has a pathwise redefinition on the space C^'^'^ 

Proof. Step 1. This prehminary step is devoted to a few details used below. Recall 
that the classical Sobolev space (M'^) is defined as the space of all f £ (W^^ having 
distributional derivative Df G (M'^^''). Recall also (see remark 4 of section 2.5.1 of 
[21] ) that, for every e E (0, 1), the space W^^^ (M°') of remark [34l is characterized as the 
space of all / e (M'^) such that 



/ 



\x — y\ 



and as a norm on W^^^ CMJ^^ one can take the following one: 



= ii/iii + mwi + 1 '-;z::r ^^^y- 



\D f(x)-D f{y)f 
\x - y\ 



Then it is easy to verify that for every e, e' G (0, 1) with e > e' the following assertion is 
true, where B (0, R) denotes the ball of center and radius R > 0: 

f G C^+^ (r'^^ , / with support in B {0, R) ^ f e W^+'' 

and 

\\f\\l^^,<C{R,e,e',p,d)\\fr^,^, (32) 

where C {R,e,e' ,p,d) is a constant depending only on R,e,e' ,p,d. 
Indeed, we have 

k-yr+^^ JB{o,R)xBio,R) |j;-yr+^P 

J \x-y\<l,\x\<R,\y\<R \x - y\ 

x-y\>Ux\<R,\y\<R \x-y\ 



< 



p 



where C {p, d) is a constant depending only on p, d. The claim (|32p easily follows from 
this inequality. 
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Step 2. Let d and e be given, as in the claim of the theorem. Choose e' E (0, e) and 
p> jr- Let W be defined on a complete probability space {Vl,!F,P). Remark [Ml states 
that there exists a random variable C > such that, for every ip S Wp~^^' (M'^), 



T 



{p{Wt),dWt: 



on a full probability set fi^. 

For every R > 0, let 0^ : ^ [0, oo) be a function such that 9ji (x) = 1 for 
\x\ < R + 1, 9r{x) = for \x\ > R + 2. Given v9 G C^+^ (M'') , we have 99 • G C^+'^ 
and thus (f-ORe W^^"' (M'^) . Therefore 



{{^■Ob) {Wt),dWt: 



<c\\^.en\\^. 



on a full probability set 0,^.0^. 

From step 1, there exists a random variable Cr > 0, independent of ip, such that 



{iip-eR)iWt),dWt: 



on ^l^.e^, 



where we have also used the fact that \\p ■ 6r\\(ji+^ < Cq \\(p\\(^x+s for some constant C^i > 
depending on the function 9 (and thus on R again). Redefine, if necessary, Cr in such a 
way that R 1— > Cr is non decreasing, with probability one. 
Let Ar be the set 

Ar = {tr > T} , 

where tr is defined in remark [33l The sets Ar increase with R. Given the family of 
events Ar and random variables Cr, we can define a new random variable C" > such 
that Cr < C on Ar (it is sufficient to put C = Cn+i on j4Ar+i\j4Ar). Thus, given 
G C^+' (R'^), we have 



< C' \\ip\\ci+e on n^.g^^n Ar. 



For every R > and (/? G C^'^'^ (M*^) , there is a P-null set Nr^^p such that 



((v5-^K)(W^t),<iVFt; 



r 



((^(W^t),dTyi) on Afl\iVfi,^. 



Therefore, given i? > and G Ci+^ (M'^) , we have 



{^iWt),dWt] 



< C \\<p\\ci+e on VL^Sr n ^i?\iVij_^. 



It follows that 



Jo 



{p{Wt),dWt, 



<C'M\ci+e on y (J^^.e^ n ^ij\iVij,^ 

fi>0 
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Since P (Uii>o Ar) = 1 we have P (U/j>o n Ar\Nr^^)) = 1. This means (Wt) , dWt) 

has a pathwise redefinition on the space C^^'^ {W^^ . The proof is complete. ■ 

Remark 36 The strategy of step 2 in the previous proof can he used to deal with function 
spaces of Frechet type that are not Banach spaces: by localization of the stochastic process, 
one can restrict the attention to compact support test functions and then prove the existence 
of a pathwise redefinition in topologies without decay at infinity. For this reason, even the 
uniformity in x in the definition of C^^^ (W^) is not necessary. 



Remark 37 In rough path theory (see U5^ ). for every rough path j of a certain class which 
includes a.e. path of Brownian motion, a notion of integral Jq {f {^t) , djt) is defined for 
every function ip with e -Holder first derivative (for arbitrary e > 0). The previous theorem 
is conceptually similar; a closer comparison, however, requires further investigation. 



Finally, we have to prove lemma 



4.1 Proof of lemma 1291 

If max[o,T] \Wt\ < |x| /2 then, for every t £ [0, T], 



and thus 



1 , , , , 3 , , 

- \x\ < \x — WA < - \x\ 

exp (— e \x — Wt\) < exp {—e \x\ /2) 

\2d-2a 



1 



\x-Wt 
1 



2d-2a 



< ^^^^j. . if 2d - 2a < 



\X - Wt\ 

exp (—£ \x — Wt 



< 



\x\ 

22d-2a 



\X\ 



\2d-2a 



if 2ci - 2a > 



X — WA \x\ 



for a suitable constant Cn d > 0. Therefore 



E 



exp {—e \x — Wt\ 



p'/2 



lo \x - Wt 



\2d-2a 



-dt 



dx<Ii + h 



[ e.p{-ep'\x\/A)E i f 

JK'i \J0 \x\ 



p'/2 



dt 



dx 



1-2 := 



E 



max[o,T] lWt|>|a;|/2 




p'/2 



\x - Wt 



\2d-2a 



dt 



dx. 
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Obviously Ii < oo, being {d — a + 1) p' < d. Moreover 



s 

l2< I p(max\Wt\>\x\/2]^^' E 

V V[0,T] ' " ' 



1 



\ {l+<5)p'/2' 

\x - Wt\ J 



dx 



for every 6 > 0. Recall the exponential inequality (see [l8] Proposition 11 



P max Wf > B ] < e 2t . 



It easily implies, by symmetry, that 

P { max \Wt\ > \x\ /2 ) < 2e~3T 

W\m J 

and thus there exist Cs^Xs > (depending also on T) such that 

s 

P ( max \Wt\ > \x\ /2 ) < C^e^^^l^L 
Moreover, by Young inequality, 

1 1 
Q-^+s = a'^+s ■ 1 < a + Cs 

for some constant Cg > 0. Thus, for every 6 > 0, 

f-T 



l2<C's + C's I e-^^l^'l^ 



1 



{l+5)p'/2- 



\x -Wt 



l2d-2a 



dt 



dx 



for some constant > 0. 

The following lemma is inspired from the proof of Corollary 2.4 of Elworthy, Li, Yor 
[2] and in fact it was suggested to us by K. D. Elworthy. 

Lemma 38 For every (i>2, g>l,0GM, xG M'*, we have 



E 



dt 



\x + Wt 



2(1-6) 









Ix + Wrf" 


+ \xf^ 


+ [ E 






Jo 



dt 



Proof. Consider the process Zt = \x + Wt\'^ (squared Bessel process of dimension d). 
From Ito formula we have 

dZt = 2{x + Wt, dWt) + dt, Zo = |xp . 
Introducing an auxiliary one-dimensional Brownian motion (/3j) we may also write 

dZt = 2y^tdl3t + dt. 
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Since d > 2, the one-point sets are polar sets for a d-dimensional Brownian motion, see 
Proposition 2.7, p. 191 of [18]. Therefore P {Zt > 0,t £ [0, T]} = 1 and we can develop 
using Ito formula for any G M. We obtain 



d[Z^ 



e/2 



2 * 



(e-2)/2 



iee-2 



2^/ZtdPt + dtj + --^Zf ^^^hZtdt 

= ezf-^^^^d(i,^cezf-^^i^dt 



where cq 



e(g-i) 



. Therefore 



QZ, 



(e-l)/2 



dUt = 4^' - zi'^ 



ceZ- 



(e-2)/2 



dt 



and thus, from BDG inequality, for every q> 1 



E 



f-T \ 9/21 \ / fT \ 

J 9'zt'dt] <c,E U ezl'-'^/'d(3A 



< Cg^e <^E Zl 
This implies, by Holder inequality. 



eq/2 



+ E 



Z\ 



eq/2 



+ E 







zf-^^"dt 



E 



1 .V 

w-'l 



z, 



eq/2 



+ E 



Z\ 



eq/2 



+ E 



1 







{2-e)q/2 



dt 



and the proof is complete. 



We go on with the proof of lemma [29l Simply by taking 1 — 9 = d—a and q = {I + S)p' 
we have: 



E 



1 



\x — Wt\ 

< cjs 

+ C I E 



2d-2a 



dt 



{l+S)p'/2 



1 







1 



+ 



|{d-o-l)(l+(5)p' 



dt. 



With the notation pt (y) 
we have 



^=^exp -f^ ,, 



E 



1 



\x + Wt\ 



(d-a-l){l+S)p' 



and the bound Pt {u) dt < '^''°^^^L2^^^^ 



1 



\x + y 



{d-a-l){ 
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and 



T 

E 







1 



1 



dt 

■T 



I , |(d-Q+l){l+5)p' , 



Pt (y) dt dy 



Thus, with a new constant C > depending on 8 and the other parameters, 

/2<C(l + /«+/f +4=^) 



where 



F + 

|(d-o-l)(l+<5)p' ' 



X 



(3) ._/■/" r-^M 1 ^^^""P^'l^'^ rfrA; 



Choose 5 > such that (d — a + 1) (1 + (5) p' < d. Since 

{d-a-l){l + 6)p <{d-a + l){l + 6)p <d 

the term 12'^^ is finite. For l!^^ and it is sufficient to integrate first in x, bound the 
result uniformly in y, then integrate in y; one proves that if^ and if ^ are finite. The 
proof is complete. 

5 The energy of a random vortex filament 

In [H [5l [7] with the purpose of modeling turbulence in 3d fiuids, the authors introduce 
and study a model of random vortex filaments based on Brownian motion. This model 
has been extended to the fBm with H > 1/2 by [16j and [8j. Here we recall briefly the 
model, emphasize the relationship of the vortex energy with the pathwise regularity of the 
current associated with the vortex core and obtain new conditions for the integrability of 
the vortex energy for the case H G (1/4, 1/2). 

For simplicity we consider only a single vortex since extension to a linear superposition 
of different vortexes is straightforward (and even to a random field of Poissonian vortexes, 
see for example f6J). Let {Xt)te[o,T]j T > be a 3d fBm with Hurst parameter H £ (1/4, 1) 
and consider the associated vector current, formally written as 



eo(x)= r 5{x-Xt)dXt 
Jo 
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where the integral is a symmetric (Stratonovich) integral. This object should be under- 
stood according to theorem [11] that is as a random distribution in the Sobolev space 
H~°'[M.'^) of sufficiently large negative order. The vorticity field is then built by superpos- 
ing translates of this core weighted according to a compactly supported signed measure p 
with finite mass which determines the intensity of vorticity. For more details about those 
considerations, the reader can consult [4]. Then we end up with 

Cix) = / Co{x - y)p{dy) 

which is again a random distribution. 

The velocity field u is generated from ^ according to the Biot-Savart relation 

u{x) = I /C(x - y) A Ciy)dy = [ /C * p(x - y) A Uy)dy (33) 

where A is the vector product in and the vector kernel IC{x) is defined as /C(x) : = 
{ATr)~^x/\x\^ and /C * p denote the convolution {]C*p){x) = ]C{x — z)p{dz). The kinetic 
energy of the fluid is then defined as the L^(]R^) norm of u: 



£ = / \u{x)\^dx = \\u\\\ (34) 

It is then interesting to find conditions on p such that the kinetic energy of the fluid is 
finite. Abstractly we have u = where we introduced an operator <I> whose kernel is 
IC* p having Fourier transform 

HIC*p){q) = ^m 

\q\ 

where we denoted p the Fourier transform of the measure p. 

From now on will stay for L^(M^). Since, by Corollary 1221 Co belongs a.s. to the 
space i?'^"(M'^) for any a > an = l/(2i^) + 1/2 (since d = 3), the condition u e L"^ a.s. 
can be satisfied if ^ : ff2^"(M^) which in Fourier variables is sufficient to require 

that ^ 

||$||^_._^, = ||$(l - A)"/2||^2_^2 = ess sup ^^(1 + |<7p)"/' < 00. 

^ geM3 1^1 

for some a > an- 

We can now formulate the following result. 

Corollary 39 The kinetic energy of the vortex filament ^ built upon a 3d fractional Brow- 
nian motion of Hurst index H > 1/4 is a.s. finite and in if the measure p satisfies 

esssup|^(g)||gr^(l + |gp)"/^ < 00 (35) 

for some a > an- 
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Remark 40 1. Known conditions on p which guarantee the integrability of the energy 
are given in for the case of Brownian motion and ltd Brownian processes, and 
in ]lb^ for the case of fractional Brownian motions with Hurst parameter H > 1/2. 
From \Td^ it can he deduced that a sufficient condition for the integrability of the 
energy is 



j 



or, written in a different but equivalent form, 

p{dx)p{dy) 



\x — y\^l^^^ 



< oo. 



2. Condition [35\) implies Condition ^36\) when H > 1/2. 
In fact the left-hand side of Condition 1135]) gives 



I dq\p{q)\W-\l + \q\''r + kl')""kl'/^"^ = / dq\p{q)\^\q\^l''~^{l + \q\ 

r+oo 

< A / r^/^n+r^ydr, 



where A is the finite quantity of [35]). Clearly previous expression is bounded for 
a > an- 

Of course the converse is not true. 

3. A way of finding similar conditions to 136\) is to follow the steps of Sec. for the 
Hilbert space regularity of the stochastic currents and rewrite (formally) the kinetic 
energy as 

£ = [ [ {dXt,g{Xt-Xs)dXs) (37) 
Jo Jo 

where g is a vector kernel with the following Fourier transform 
and Hq is the following matrix 

(^q)al3 = ^al3 — y"^; a, (5 = 1, ... ,3, 



which projects in directions orthogonal to q. Formula [31) can be understood, for- 
mally according to Theorem \ll\ as being the limit of the expectations of e- approximations 

£e = r r 9{Xt - Xs) {DeXt, D,X,) dtds 
Jo Jo 

To obtain conditions for its finiteness in the spirit of eq. 136]). we need to follow again 
the computations involved in the proof of Theorem{T^ and use a different strategy in 
bounding some terms. Then we can prove the following: 
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Theorem 41 Let H > 1/4 and let p : M.^ ^ be a function with Fourier transform p 
satisfying 

then the family of random fields {ue}s<^(o,i) defined as 

Ue[x)= [ {IC*p){x-Xt)AD^Xtdt, xeM.^ 
Jo 



converges a.s. in ^(17; L^(M'^;M^)) for any 9 > to a random field u £ L 

Proof. We will prove that sup^g^g i) -EUn^lp < co following the lines of the proof of 
Theorem 1181 then the conclusion follows applying Theorem [TTJ Let = ||'Ue|P- 
We start by treating the case H > 1/2. 

Using Theorem [23] (Wick theorem) and independence of different coordinates we have 

E£e = E r CV] gu {Xt - X,) CoviD^.X^ D^.X}) dtds 
Jo Jo ^ 

r TYI V^V,5i,• {Xt - X,) \Cow{Dlx},Xj - dtds (39) 

= E f [ Trg {Xt - Xs) Cov(D°X/, L>°X]) dtds 
Jo Jo 

since a direct computation shows that J2i ^igik{x) = 0. Using the first bound in Lemma [T7l 
we get 

E£e<constE [ [ Trg {Xt - Xs) \t - s\'^"~'^ dsdt 
Jo Jo 

= const / dqTrg{q) [ [ \t - s\'^^-'^Ee-'^'J'^'-^^^ dsdt 
Jr^ Jo Jo 



Jo 



= const / dqTTg{q) / / 
Jr^ Jo Jo 

< const / dq\g{q)\ / dt / drr^'^-'e-^''^"-"" 
Jr^ Jo Jo 



(40) 



const / dq\g{q)\\q\^/^-^ dt dyy^-^^^e'y^ 

Jo Jo 



T foo 



< COnstT / dq\g{q)\\q\^/^-^ 
JR3 



since 



oo 

2 , 
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for 1 — 1/H > —1, that is H > 1/2. Sufficient condition for uniform boundedness of ES^ 
is that 

' dq\g{q)\\q\'/''~'^ < 



oo. 



Let us now consider the case H < 1/2 and rewrite the approximated energy as 
"1 rl rl rl 



= \\Uef 



[ [ {h{Xt-Xs)DeXs,D,Xt)+ [ [ {g{0)DeXs,DeXt 

Jo JO Jo Jo 



where h{x) = g{0) — g{x) > 0. Note that g{0) is weh defined using the hypothesis of the 
theorem about the integrabihty of its Fourier transform, moreover, as in Thm. [T2] (in the 
H < 1/2 part) we have the hmit 

/' f\g{0)D,Xs,DeXt) ^ {{Xi-Xo)g{0),{Xi-Xo)). 
Jo Jo 

So let us focus on the double integral with the kernel h. Proceeding as in the H > 1/2 
case we have 



J = - [ [ {h{Xt - Xs)D,Xs, DeXt) 
Jo Jo 

< -const / f [ dqTi:h{q)E[e'^'i'^'-^^^\t - s\'^"-^ dtds 
Jo Jo Jr^ 

= -const r r [ dgm(g)e-l^l'/2(*-)'''|t-s|2^-2citcis 
Jo Jo Jr3 

but since h{q) = g{0)S{q) — g{q) we have 
/ dgTr/i(g)e-l«l'/2{t-^)'^ = [ dqTr[g{0)6{q) - 5(g)]e^l'?l'/2(t-^)'" 

dqTr[g{0)5{q) - ^(g)] [g-l^l'/^C*"^)^" - 1] = - / dqTTgiq)[e-\'^\' 



Then 



|J| < const r r [ dq\g{q)\il - e-l^l'/2(*-^)''')|t - sl^""'^ dtds 
Jo Jo Jr^ 

< const/ dq\g{q)\ I dt [' dsit - sf^-\l - 6-^^"^'^'-''^'" ) 

Jo Jo 

< const / dq\g{q)\ dt drr^^'-^il - e-^'^^'^'" 

Jr3 Jo Jo 

= const / dq\g{q)\\q\'/'''' / dt / dyy'-'/''{l - e'^'l^ - 1) 
Jm? Jo Jo 

< const / 

Jr? 
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where we made a change of variables y = \q\T^ and we used the fact that 

/ dyy^-^'^H - e-y^'^) < / dyy^-^'^ min(y2, 1) < oo 
Jo 

since —3 < 1 — \/H < —1 when 1/4 < H < 1/2. So we obtain the uniform boundedness 
of E£^ when eq. ([38|) is satisfied. 

Analogously, the case H = 1/2 does not pose any additional problem. ■ 

Remark 42 Note that for H > 1/2 we recover condition i36\) . However while in Q {Wjj 
only the existence and the integrahility properties of the energy are studied, here we have 
also informations about convergence of e- approximations of the velocity field generated by 
the random vortexes. 

A Some proofs and auxiliary results 

Proof of Lemma [2l Denote , ^/ ^^/g by 7, for shortness. Notice that for x = we 

r(a)(47r) ' 

have 

Ka (0) = 7 / t"~ie"*— < 00 if and only if a > -. 

Jo t 2 

2 

For X 7^ we may use the change of variables t = |x| s and get 

Ka{x) = Ixl^""*^ /> (x) 



/9 (xj := 7 / s 2 e 4s 

Jo ^ 

where the integral converges for every value of the parameters, thanks to the exponentials. 
For < a < |, we have 

< 7 /"%"-fe-i-l-l'^- < p(x) 
Ji s 

for a positive constant Ca^d- Moreover, 

/ X / a-i -±ds f°° Q,_d _W_u|2 (is 

/9 (x) < 7 / s 2e 4s s 2e i i^^i — 

Jo s Ji s 

\x\ 

< Cn.d I e »s— + Ca,de 4 / s° 2 — 

. J_ s 







<Ca,de +Ca,de '4 |x|2 "<Ca,de 's' 
for a positive constant Cq d that we do not rename at every step. 
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If a > ^, we directly have from the original formula 

c'„,rfe-— < 7y r-ie-— < (x) 

<7/ t 2e 4t 1-76 4 / t"2e'^ — 

Jo * J|x| ^ 

< c;,e-^ 

for a positive constants c'^ ^, ^. To estimate Kq.(0) — A'q(^) we write 

i^„(0) - K^{x) = 7 / e-*(l - e-^)^ 

JO * 

= 7|xr-i rs"-fe-^N'(l-e-i^)- 
Jo s 

and use the same arguments as above. In the case a = d/2 we simply split the integral 
as above and by straightforward estimation we can prove that Ka{x) < const log|x| for 
small |x| and that Ka{x) decay exponentially for large The proof is complete. ■ 

Proof of Lemma [4l We observe that, since Ka is the kernel of the operator (1 — A)~" 
we have the identity Ka~i{x) = (1 — /S.)Ka{x) so that 

-/^Ka{x) = Ka-l{x) - Ka{x) XJ^O. 

Then | — AKa{x)\ < |-K'q,_i(x)| + which gives the required bound using Lemma [2l 

■ 

Proof of Lemma 1171 We start with the first estimate in 1. A direct computation shows 
where 

|l + ^|2// + |i_a.|2H_2 
^(^) = -2 • 

The function <I> is continuous in (0,oo), lim^.^o '^'(2^) = 2ff — 1 so, when |t — s| < 2e 
we have ^ 

~ $ I 1 < const \t - 



2 \t-s 

Moreover lim^^^ioo \x\'^~'^^ (^{x) = 2, so when |i — s| > 2e there exists a constant not 
depending on e such that 

~ -^> f < const £^^-2 < const It - sP-^"^ 



2 \t- s 

which proves the first claim. 
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We discuss now the second estimate in 1. and point 2. A direct computation gives 
CoviD«,XlXl-Xl) = -Cov{D»X',.Xl-Xl) = ^^{\t-s + 4'« -\t-s-E\'«) 

where ^^{x) = ^-^-^ — 2x^~^^ — • I* is easy to show that ip is continuous and ■0(0+) = 
moreover \x\'^~'^^ ip{x) — > 2H when x — > ibcxD. This allows to conclude the proof. ■ 

Proof of Lemma I27L The first estimate in 1. is very similar to the previous lemma. A 
direct computation shows 

Cov{D;XI,D-XI) = I ^ $ (^__ j (41) 

where is the same as previously. As for the other points a direct computation gives 
Cov(D-X;, X; - XI) = -Cov(D-Xl,Xl - Xj) = i (|t - « + e\™ - |( - s - E\™ - E™) 

- «-.r-'*(^) 

where ■(/'(x) = - — ^^~2x~^^ — • Then, it is not difficult, arguing as in lemma [T71 to conclude. 
In particular one can evaluate the limit in (I4ip . 
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